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Quantum simulations based on near-resonance Bose gases are limited by their short lifetimes 
due to severe atom losses. In addition to this, the recently predicted thermodynamical instability 
adds another constraint on accessing the resonant Bose gases. In this Letter, we offer a potential 
solution by proposing long-lived resonant Bose gases in both two and three dimensions, where the 
conventional few-body losses are strongly suppressed. We show that the thermodynamical properties 
as well as the lifetimes of these strongly interacting systems are universal, and independent of short- 
range physics. 


Cold atomic gases near resonance offer an ideal plat¬ 
form for the study of strongly interacting quantum sys¬ 
tems and have recently attracted considerable interest. 
It not only provides a distinctive approach to advanc¬ 
ing our understanding of many puzzling quantum many- 
body problems in various other fields of physics, but also 
offers intriguing topics of its own. One outstanding ex¬ 
ample is the atomic Bose gas near resonance, which is 
experimentally accessible thanks to the extreme tunabil- 
ity of interatomic interactions. This allows for experi¬ 
mental studies on Bose gases at large positive scattering 
lengths Alongside these experimental works, vari¬ 

ous numerical and theoretical efforts have also been made 
to analyze this strongly interacting system Ifll- 221. 


One of the main constraints that limits the experimen¬ 
tal access of resonant Bose gases is their short lifetimes 
due to few-body losses. These are recombination pro¬ 
cesses which lead to the formation of few-body bound 
states and consequently atom losses. The minimal num¬ 
ber of atoms involved in these processes is three, and the 
corresponding process is known as three-body recombi¬ 
nation. There have been extensive studies of this pro¬ 
cess in the few-body context [1, 0, ■ It has been 

shown that, apart from interference effects, the three- 
body recombination rate tends to increase rapidly when 
resonance is approached, which significantly reduces the 
lifetime of the gas. This is an intrinsic issue of atomic 
Bose gases which practically limits the outcome one can 
get out of this platform near resonance. Overcoming this 
difficulty has become one of the main focuses of recent 
experimental studies. For instance, a very interesting 
attempt applied quenching of the scattering length to 
access resonance and a local equilibrium was reached Q . 


Another complexity of resonant Bose gases is the re¬ 
cently predicted thermodynamical instability near reso¬ 
nance |l2l . [iM3, 111- Different from few-body losses, 
this issue of thermodynamical instability is a many-body 
effect. It is characterized by a sign change of the com¬ 
pressibility at a critical scattering length, beyond which 
the compressibility becomes negative. In two dimensions 
(2D), a diagrammatic analysis further implies an anoma¬ 
lous behavior of the compressibility as a precursor of the 
instability 0, which is consistent with a recent exper¬ 


iment on 2D Bose gases [^. In three dimensions (3D), 
however, there is still lack of conclusive experimental ev¬ 
idence for this instability. Nevertheless, with this pre¬ 
dicted many-body instability, accessing resonant gases 
might run into a sharp phase transition rather than the 
smooth variation of lifetime predicted by the few-body 
losses 0. 

Despite these difficulties, efforts have been made to ex¬ 
plore the universal aspects of the gas iii. Naively, 
the only relevant length scale of the system at resonance 
is the interatomic distance, which should universally dic¬ 
tate the low-energy physics of the system. However, for 
resonantly interacting bosons, there are a series of infinite 
three-body bound states with a discrete scaling symme¬ 
try, i.e., the Efimov states 0. The binding energies of 
the Efimov states further depend on the short-range de¬ 
tails of the interatomic potentials. Therefore, it has been 
proposed that the many-body properties of the gas would 
further depend on the non-universal features of the in¬ 
teraction 33“32|. From the few-body point of view, this 
underlying Efimov physics manifests itself in the three- 
body recombination process as the interference peaks and 
dips in the three-body loss rate [0. This has been uti¬ 
lized to detect the Efimov effect in experiments [ii-ili. 
On the other hand, in the many-body limit, so far there 
has been no indisputable experimental evidence of the 
Efimov effect in resonant Bose gases 0, @| . More exper¬ 
imental efforts might be needed to address the issue of 
universality. 


Given the current challenges of studying resonant Bose 
gases, it is tempting to ask whether there exist long-lived 
resonant Bose gases where few-body losses are strongly 
suppressed and that are also thermodynamically stable 
and universal. In this Letter, we propose such a sys¬ 
tem which can be potentially realized in optical lattices 
with long-range hopping. It provides a stable degenerate 
quantum gas which is an ideal candidate for quantum 
simulations of the thermodynamics and the universali¬ 
ties of strongly interacting systems. 


We hrst consider a 3D Bose gas with a power-law 
single-particle dispersion e(fc) ^ {S is a small 

number) and contact interatomic interactions. It can be 
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described by the grand-canonical Hamiltonian 
H = ^(e(fe) - n)b\hk + 2Uono ^ blbk 


UqUo ^ bkb-k 


k,q 




q + h.C. 


+ ^ H +^-C- (1) 


2H 


.q 


where (6^) is the bosonic creation (annihilation) oper¬ 
ator, with no the number density of the condensed atoms, 
/i the chemical potential of non-condensed atoms, and 
Uq the strength of the contact interaction. H is the vol¬ 
ume of the system and the summations are over non¬ 
zero momentums. We set the reduced Planck constant 
to be unity and choose units so that e(fe) = fc^/^+^/2 
where k = |fc|. This kind of dispersion can be realized 
in optical lattices with long-range hopping. Here we take 
the 3D cubic lattice as an example and consider a long- 
range hopping varying as Uj ^ \ri — , where 

tij denotes the hopping amplitude between lattice sites 
i and j. This long-range hopping leads to a dispersion 
e(fe) ^ near the bottom of the lowest band. There¬ 

fore, the system described by Eq. © can be realized by 
loading bosons into such a lattice, with the atom number 
per lattice site being much less than one. The tuning of 
interactions can be achieved by a resonance driven by the 
effective mass, a standard method for keeping the system 
in the single-band regime [s^. 

One direct consequence of a general power-law disper¬ 
sion is the breaking of Galilean invariance, which con¬ 
siderably alters the two-body scattering behavior. The 
T-matrix for a scattering process with total momentum 
k and total energy E takes the form 


T{E,k) 


4^2j[l+0(d2)] 

- {-E - i0+)“ -f f{E, k )6 ’ 


( 2 ) 


where B is the binding energy of the two-body bound 
state, a = |y|^, and / is a regular function satisfy¬ 
ing f{E,0) = 0. One can immediately see that T{E,k) 
possesses a completely different form compared with the 
Galilean invariant case, where T{E,k) is simply given 
by T[E — fc^/4,0). Furthermore, the fe-dependence of 
T{E,k) here is suppressed by an extra power of S and 

is thus irrelevant when (5 <C 1. It is convenient to define 

_ 1 

an effective scattering length X = B 3 / 2 +s ^ which dic¬ 
tates the low-energy scattering. It renormalizes the bare 
parameter Uq via 


r (Pk 1 

J (27r)3 k3/2+s 


1 

4^2y3/2-5j[l-^0(J2)]' 


(3) 


The divergent A marks the position of resonance, which 
is the threshold to form two-body bound states. 


To analyze this interacting Bose gas, we apply a self- 
consistent approach, where the energy density of the sys¬ 
tem, E(nQ,fj,), for a given uq and /r is first calculated. 
The chemical potential of a gas with atom density n 
can then be calculated using the following set of self- 
consistent equations. 


dE{no,p 

k-cijiQ, p = ---, = fJ,c{no, p 


n = no 


dno 

dReE{no, fj.) 


(4) 


where fXc is the chemical potential of the condensed 
atoms, which should be equal to p, in equilibrium. The 
second line of Eq. (|4]) is the number equation. The cal¬ 
culation of E(no,/i) is carried out diagrammatically us¬ 
ing the effective field theory method [3, [3 IHj. It 
contains contributions from all A^-body diagrams with 
N = 2, 3, • • •. In the diagrammatic calculation, instead of 
Uq, the T-matrix parametrized by A (Eq. ((^l) is adopted. 

We first discuss the weakly interacting or the dilute 
limit, which is defined by the condition X^nS 1. The 
dimensionless parameter X^nS is an analogue of the gas 
parameter Pn for 3D Bose gases in free space, where a 
is the scattering length. In the dilute limit, the lead¬ 
ing order contribution to /i is given by the Hartree-Fock 
chemical potential /thf = 47r^A^/^“'^nd. The dilute con¬ 
dition states that /thf is much smaller than B, which 
makes it possible to carry out an expansion in the di¬ 
lute limit. More specifically, the next-to-leading-order 
contribution to the chemical potential is proportional 
to /rHF(A^nd)“. It comes from the ^-dependence of T, 
which results from the energy shift of the non-condensed 
atoms due to many body effects. For the quadratic dis¬ 
persion with 5 = 1/2, this correction is proportional to 
/^hfVA^ n, resemblin g th e Lee-Huang-Yang correction for 
dilute Bose gases [3J, . The next order correction con¬ 

tains an extra factor of 5 but gives the leading order con¬ 
tribution to the imaginary part of the chemical potential, 
which scales as p,HF(A^ni5)i5. This imaginary part comes 
from the pole structure in the T-matrix and reflects the 
atom losses due to three-body recombination processes. 
Gompared with the real part of the chemical potential, 
the imaginary part is further suppressed by a factor of d, 
in addition to the expected suppression by the diluteness 
parameter. This indicates an extremely low few-body 
loss rate and long lifetime. 

By solving the self-consistent equations using the 
method of iteration, we obtain the solution: 


M = Mhf 


,3 81n2 — 4 1 2 ^3 

1 + ( - +-^-d ) {SPX^nSr 


—iA'K{d>PPn5)6 -I- 0((X^nS)S^) 
1 


no = n 


1 - ^(8PX^nS)°^ + 0((X^nS)S) 


(5) 
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The expansion in Eq. breaks down near resonance 
when A^n(5 ~ 1 or ^ 1. However, the presence of the 
small parameter S makes it possible to systematically an¬ 
alyze interacting Bose gases even at resonance. Follow¬ 
ing similar techniques developed in Ref. 21|, we carry 
out this analysis by examining all the diagrams that con¬ 
tribute to the energy density (Readers who are not in¬ 
terested in these technical details can skip to Eq. ([9|)). 
Each diagram can be described by two integers, N and 
L, which are half the number of external condensed 
atom lines and the number of closed loops respectively. 
The corresponding diagram is called an iV-body L-loop 
diagram. When X^n6 ^ 1 or ^ 1, the energy de¬ 
pendence of the T-matrix begins to play an important 
role. This makes the self-consistency necessary even at 
the tree level. In the tree level, the energy density is 
= nQT{2fi, 0)/2, where the superscript of E 
indicates the number of loops. Consequently, the self- 
consistent equation reads. 




47r^no(5[l -|- 0(d^)] 


A-(3/2-5) _ (_2^_i0+)" 


( 6 ) 


In the limit when A^ni5 ^ I, the solution of Eq. takes 
a simple form. 


= (27r2no(5)^/^+'^/^[l + 0(<5)]. 


(7) 


Beyond the tree level, one can easily show that all dia¬ 
grams are free of ultraviolet divergences after the energy 
dependence of the T-matrix is taken into account. To an¬ 
alyze the contributions of these diagrams, we take as 
the leading order result for the chemical potential. The 
contribution to the chemical potential from an 7V-body 
T-loop diagram can then be estimated to be 

( 8 ) 

which only depends on L, as indicated by the superscript 
of /r on the left hand side. One can see that, compared 
with the tree-level result, contributions from diagrams 
with T > 0 are suppressed by extra factors of S, which 
validates taking as the leading order result. More¬ 
over, Eq. ® shows that the contributions from diagrams 
beyond tree level are organized according to the number 
of loops, independent of N. This suggests a systematic 
expansion in terms of L controlled by 6 , even though the 
gas is near resonance or X^nS 1. 

In the limit X^n6 —>■ oo, we solve the self-consistent 
equations rigorously up to order 0{S) by summing up 
all the one-loop diagrams. The self-consistent equations 
lead to a solution that is universal, 

2v2 3/2+5 o 

^ = -^eF5^^(I + 1.225-^3.795-1-0(5^)) 
no =-n.(l-(-0.405-I-0(5^)) (9) 

O 

where ep = (67r^n)^/^“''^/^/2 is the Fermi energy defined 
for the Bose gas with the single-particle dispersion e{k) = 


TABLE I. The chemical potentials of Bose gases with disper¬ 
sion e(fe) = A:‘*0+^y2 (d is the spatial dimension and 5 is a 
small number) in the dilute limit and near resonance in both 
two {d = 2) and three dimensions (d = 3). Note that they 
have the same form in both 2D and 3D. The dimensionality- 
dependent numerical coefficients, Di, D 2 , D 3 , Ri, R 2 , and 
Rs, are shown in Table IIIl Here (3 = , /thf ~ X'^^^~^n5 

represents the Hartree-Fock chemical potential, and -q = X'^nS 
is the dilute gas parameter, with A the effective scattering 
length and n the atom density, tp ~ jg the Fermi 

energy defined for density n. 



Rep 

Imp 

Rep 

Dilute limit 

in < 1) 

Phf[1 -f Diq^ -1- D2q^5 H- ] 

—DsqS + ■ ■ ■ 

Resonance 

d/2+5 


{q 00) 

epl-tiid ‘‘ (i-I--f • ■ • )1 

—R3S - 


TABLE 11. Numerical coefficients in solutions of the chemical 
potential in 2D and 3D. The solutions with these coefficients 
are shown in Table HI 



2D 

3D 

Di 

127r 

127r2 

D 2 

—87r(2 -I- 5 In 2 -I- 3 In tt) 

16f g +21nF) 

D 3 

48tt^ 

32-k^ 

Ri 

273 

272 

3 

3 

R 2 

2.08 

1.22 

Rs 

5.89 

3.79 


^3/2-i- 5^2. It can be seen that the imaginary part of the 
chemical potential is smaller by a factor of 5 compared 
with its real part, implying a long-lived Bose gas near 
resonance. We shall discuss this in more detail below. 

We now carry out a similar analysis for a 2D Bose gas 
with a nearly-linear dispersion e(fe) ^ , which yields 

similar results. In the dilute limit, by defining a 2D gas 
parameter A^n5, the expansion of in terms of the 

gas parameter and 5 has the same form as Eq. ® except 
for different numerical coefficients. At resonance, this 
2D Bose gas is also shown to be universal and long-lived. 
Furthermore, apart from different numerical prefactors, 
the result expressed in terms of the Fermi energy is orga¬ 
nized exactly in the same way as Eq. ®. Table U summa¬ 
rizes the main results for both cases and the numerical 
coefficients are shown in Table IE 

It can be seen that the equations of state at reso¬ 
nance have the same universal form in both 2D and 3D. 
Since resonance corresponds to a fixed point of the run¬ 
ning two-body coupling constant, this implies that the 
3D Bose gas with dispersion e(fc) ^ and the 2D 

Bose gas with dispersion e{k) ^ k^~^^ belong to the same 
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universality class near resonance. This can be further 
generalized to a series of Bose gases with e(fe) ~ 
where d is the dimensionality. A detailed discussion of 
this universality class in different spatial dimensions will 
be presented in an upcoming article. 

One important feature of the systems discussed above 
is their long lifetimes. The lifetime of the gas is deduced 
from the imaginary part of the chemical potential, which 
can be generally applied both in the dilute limit and near 
resonance. This can be understood by considering the 
condensate wave function whose time evolution 

is given by A negative imaginary part of the chem¬ 

ical potential implies a damping effect and the lifetime of 
the condensate is given by |Im^|“^. In the dilute limit, 
the leading order contribution to Im^ comes from the 
three-body scattering processes and can be directly con¬ 
nected to the three-body loss rate. The characterization 
of the lifetime using Im/i then reduces to the standard 
description in terms of the few-body losses. In contrast, 
near resonance, contributions from all A^-body processes 
with N = 3,4, 5, • • • are expected to be parametrically 
inseparable. As a result, the relation between the life¬ 
time and the three-body loss rate cannot be established 
straightforwardly. In this case, the imaginary part of the 
chemical potential is due to collective effects and natu¬ 
rally encodes the information of the lifetime. Its real part, 
on the other hand, serves as the characteristic many- 
body energy scale. Therefore, a long-lived Bose gas can 
be defined by |Im(/r)/Re(^)| <C I, which guarantees the 
realization of equilibrium before significant losses occur. 
According to Eq. ([5|) and the third column of Table. U 
|Im(/i)/Re(/i)| ^ 5, which is parametrically small. This 
not only indicates that those gases have long lifetimes at 
resonance, but also justifies our method since equilibrium 
is necessary for the self-consistent approach. 

Apart from the short lifetime, the thermodynamical 
instability that has been pointed out for 2D and 3D Bose 
gases with normal dispersions 3, 18, 281 is not present 
in the systems we proposed here. This can be seen by 
examining the compressibility, which can be calculated 
via K = {dKe^/dn)~^. Table |T] and Table HIl show that 
Rep, ^ ^/n and Rep > 0, which leads to a positive com¬ 
pressibility near resonance, indicating the thermodynam¬ 
ical stability of the gas. 

These long-lived and stable resonant Bose gases are 
universal, independent of microscopic details of atomic 
interactions. This is shown by the equation of state 
(Rep(n)) which only depends on the Fermi energy, an 
energy scale set by the interatomic distance. Moreover, 


the dynamics, i.e., the lifetime (|Imp| 


is also 


universal, again only depending on the Fermi energy or 
the density. This aspect is coincidentally consistent with 
the experimental result in Ref. Q, where the timescale 
of the dynamics of a 3D Bose gas at resonance had been 
scaled in terms of the Fermi time, tp = Therefore 
our results might shed further light on future studies of 


the dynamics of Bose gases. 

On the other hand, the systems are strongly interacting 
near resonance. This can be seen from the depletion frac¬ 
tion of the atoms which takes a universal value of I/3 and 
is not parametrically small. Additionally, the scaling be¬ 
havior of the chemical potential near resonance (p oc v^) 
is distinctly different from the linear dependence (p oc d) 
in the dilute limit where the systems are weakly interact¬ 
ing. 

In summary, we demonstrate that long-lived resonant 
Bose gases can be potentially realized in optical lattices 
with long-range hopping in both two and three dimen¬ 
sions. We show this by proposing Bose gases that have 
parametrically long lifetimes and are thermodynamically 
stable at resonance. We further prove that they are 
universal and strongly interacting, which makes them 
promising candidates for simulating universal physics of 
strongly interacting many-body systems. 

This work is in part supported by the Canadian In¬ 
stitute for Advanced Research and the Natural Sciences 
and Engineering Research Council of Canada. 
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